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ESTIMATING  .JOINTLY  SYSTEM  AM'  COMPONENT 
RELIABILITIES  USING  A  MUTUAL  CENSORSHIP  APPROACH 

by 

Hard  Doss,  Steven  Freitag,  and  Frank  Proschan 

ABSTRACT 

Let  F  denote  the  life  distribution  of  a  coherent  structure  of  independent 
components.  Suppose  that  we  have  a  sample  of  independent  systems,  each  having 
the  structure  p  .  Each  system  is  continuously  observed  until  it  fails.  For  every 
component  in  each  system,  either  a  failure  time  or  a  censoring  time  is  recorded. 

A  failure  time  is  recorded  if  the  component  fails  before  or  at  the  time  of  sys¬ 
tem  failure;  otherwise  a  censoring  time  is  recorded.  We  introduce  a  method  for 
finding  estimates  for  F  (t) ,  quantiles,  and  other  functionals  of  F,  based  on  the 
censorship  of  the  component  lives  by  system  failure.  We  present  limit  theorems 
that  enable  the  construction  of  confidence  intervals  for  large  samples. 
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1.  INTRODUCTION’  AND  SUMMARY. 

Consider  a  system  of  independent  components  labeled  1  through  m.  We  assume 

that  the  system  forms  a  coherent  structure,  which  we  denote  by  i>.  Specifically, 

the  system  and  each  component  are  in  either  a  functioning  state  or  a  failed  state, 

and  the  state  of  the  system  depends  only  on  the  states  of  the  components;  see 

Barlow  and  Proschan  (1981,  Chapters  1  and  2)  for  definitions  and  basic  facts 

relating  to  coherent  systems.  Let  F ^  denote  the  life  distribution  of  component 

i,  1=1,  2,  ...,  m,  and  F,,  or  simply  F,  denote  the  life  distribution  of  the 

<j> 

system. 

Suppose  that  we  have  a  sample  of  n  independent  systems,  each  with  the  same 
structure  i .  F.ach  system  is  continuously  observed  until  it  fails.  For  every 
component  in  each  system,  either  a  failure  time  or  a  censoring  time  is  recorded. 

A  failure  time  is  recorded  if  the  component  fails  before  or  at  the  time  of  sys¬ 
tem  failure.  A  censoring  time  is  recorded  if  the  component  is  still  functioning 
at  the  time  of  system  failure.  From  these  failure  times  and  censoring  times  we 
wish  to  estimate  F. 

In  order  to  distinguish  between  components  and  systems,  we  index  systems 
with  the  letter  i  and  components  with  the  letter  j;  i  ranges  over  1,  ...,  n,  and 
j  over  1,  . . . ,  m.  All  random  variables  are  non-negative.  We  define  the  following 
random  variables: 

T.  is  the  lifelength  of  system  i, 

X  is  the  lifelength  of  component  i  in  system  i, 

for  each  i,  X  .  .  ,  X...  ...,  X  .  are  i id  '  F . , 
lj  2)  nj  .1 


.  .  =  I (X .  .  -  T.  )  ,  where  I i A  1  is  the  indicator  function  of  the  set  A. 
l  j  li  i 


and 


Z ^ .  records  the  time  on  test  ot  component  j  of  system  i,  and  6.^  indicates 

whether  component  j  in  system  i  is  uncensored  (6.^  =  1)  or  censored  (6^  =01. 

For  each  j,  I..,  ....  Z  .  are  iid  with  distribution  H . .  The  sequence  { (Z . . ,£ . . ) ; 

1 j  ni  j  n  ij  ij 

i  <  i  s  m,  1  <  i  <  n1  contains  all  the  information  used  in  estimating  F,  and  thus  is 
called  zke  sarrZe  information. 

The  system  life  distribution  F  can  be  estimated  by  the  empirical  estimator 
pcmP  defined  for  t  5  0  bv 


FemP(t)  =1  l  KT  st). 
i=l 


(1.1) 


F  "  does  not  fully  utilize  the  sample  information.  Specifically,  it  does  not 
use  the  following  information:  the  identity  of  the  components  still  functioning 
at  system  failure  time,  and  the  failure  times  of  the  components  failing  before 
system  failure  time. 

We  propose  an  estimator  F  of  F  that  uses  the  information  described  above. 
The  estimator,  described  next,  is  based  upon  estimators  of  the  component  life 


distributions  F.,  ....  F  .  For  each  j,  let 


'(pi  :  i 


be  the  ordered 


values  of 


■lj’  “2 j 


,  Z  . .  Define 
n  i 


5(i).i  '  i 


1  if  Z  .  .  corresponds  *•.  o'  :r 

( l !  )  1 


!  0  if  Z  .  corr  spond 
(ill 


.  ;  ;  t  •  t  ,  -.e 


(When  an  uncensored  and  a  censored  r.  >’  • 


vation  is  considered  to  have  occurre 


estimator  of  F : 


F  .  1 1  >  =  1 

i 


The  definition  above  differs  from  the  usual  d>-ti: 


mator  in  that  F.ft)  is  not  arbitrari  lv  defined  t  >  *• 
J 

For  each  coherent  structure  t  of  independent 


I  III  llll  |M«  MIIM'BIP. 


a  function  h  ,  called  the  reliability  function,  such  that 


K 


pa 


•  4 


f  •  * 

?.  : 
r„v 


f  ii 


F  (t )  =  h  ( F  (t) . F  (t) )  . 

<p  <p  1  m 


(1.4) 


Here,  F>>(t)  =  l-F  (t)  and  F.(t)=l-F.(t).  A  more  detailed  description  of  relia¬ 
bility  functions  is  provided  in  Chapter  2  of  Barlow  and  Proschan  (1981).  The 
estimator  F  is  defined  by 

1  -  h.(F.  (t) ,  .  .  .  , F  (t) )  if  t<T, 


F(t)  = 


V  1 


U 


if  t  >  T 


(n) 
(n)  ’ 


(1.5) 


Here,  T^  =  max  (T^  ,T_, ,  .  .  .  ,Tn)  •  The  estimator  F  has  obvious  intuitive  appeal. 

The  properties  of  the  Kaplan-Meier  estimator  have  been  studied  extensively 
by  various  authors.  Under  the  assumption  that  the  censoring  variables  and  the 
lifelengths  are  independent,  the  Kaplan-Meier  estimator  is  the  maximum  likelihood 
estimate  (Kaplan  and  Meier,  1958;  Johansen,  1978).  Regarded  as  a  stochastic  pro¬ 
cess,  it  is  strongly  uniformly  consistent  (Fbldes,  Rejto,  and  Winter,  1980)  and 
converges  weakly  to  a  Gaussian  process  (Breslow  and  Crowley,  19~4,  Aalen,  1976, 
and  Gill,  1985)  . 

The  main  results  of  this  paper  can  now  be  stated.  Let  D[0,T1  be  the  space 
of  all  real  valued  functions  defined  on  [0,T]  that  are  right  continuous  and  have 
left  limits,  with  the  Skorohod  metric  topology.  Dm[0,T]  denotes  the  product 
metric  space. 

THEOREM  1.  Suppose  F, ,  F_,  ...,  F  are  continuous,  and  let  T  be  such  that  F .  (T)  <  1 
-  1  2  m  j 

for  i  =  1 ,  2 ,  .  .  . ,  m .  Th en  as  n  •+  ® 

K  -F0,  ....  F„-FJ-(W1(  W, . WJ 


mm  1 


weakly  in  D  [0,T],  where  W ^ ,  ....  IV  are  independent  mean  0  Gaussian  processes. 

The  covariance  structure  of  W  is  given  by 

tl  dF.(u) 

CovflV  ftp.lv  ft7))  =F  (t  )l;  ftn)  / - J— -  forOa  -t^T.  (1.6) 

■  •  “  0  n  (u)F  (u) 


F'T' 


'T  '  TV  'J  1  • 


f. 


tei 


:■  • 


r-v 

•  4 


ii  -. 


-  4  - 

Since  in  general  the  dependence  among  the  F ■ 1 s  may  be  complex,  Theorem  1  is 

not  a  trivial  extension  of  the  corresponding  result  for  the  individual  Kaplan- 

Meier  estimators  F.. 

l 

The  next  theorem  gives  a  central  result  regarding  the  estimator  F. 

THEOREM  2 .  Suppose  F ^ ,  F0,  ....  F^  are  continuous,  and  suppose  T  is  such  that 
F  j  (T)  <1  for  j  =  l,  2,  m.  Then  as  n 

n*(F-F)~*W  weakly  in  D[0,T], 

where  W  is  a  mean  0  Gaussian  process  with  covariance  structure  given  by 


Cov(W(t.),W(t,))  =  l 
j  =  l 


ah . 

*£ul- 

J 

.  ,u  ) 
m 

\ 

<U1 . um>  - 

, 

. 

3hA 

1 

•  * 

■  *  - 

^ (u  , . . . ,u  ) 
3u .  1  m 

J 

1 

(U1 . um)  = 

(fpt,) . 

(!•") 


F  (t  )F  (t0)  /  7 
J  1  J  n  u 


tl  dF.(u) 


0  Hj  (u)  F (u) 


for  OstjSt^sT. 


The  commonly  used  estimate  of  the  variance  of  the  Kaplan-Meier  estimate  is 
given  by  Greenwood's  formula  (see  Chapter  5  of  Miller,  1981).  Since  this  esti¬ 
mate  is  known  to  be  consistent  (see  Hall  and  Wellner,  1980),  it  follows  that  for 
fixed  t,  the  variance  of  F(t)  given  in  Theorem  2  can  be  consistently  estimated. 
This  enables  the  construction  of  confidence  intervals  for  Fft)  in  large  samples. 
The  method  of  estimating  F  proposed  here  has  an  additional  advantage:  the 

estimates  F  ,  F„ ,  ....  F  can  be  used  to  estimate  the  life  distribution  of  any 
1  2  m 

structure  ■]>  whose  components  form  a  subset  of  the  components  of  $  .  Specifically, 
if  h  is  the  reliahilitv  function  of  ijj,  then  the  estimate  F  defined  bv 

■j  ■  o 

F  ( t )  =  1  -  h  ,  (F.  (t)  ,  .  .  . ,P  (t) )  ,  (1.8) 

’Ll  V  l  <n 


f':'' 


■  I, .a  Mfc  ■  mim'mmu 


:i 


(1.8) 
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when  suitably  normalized,  converges  to  a  Gaussian  process;  this  fact  will  be 
clear  from  the  proof  of  Theorem  2. 

In  the  large  literature  on  point  and  interval  estimation  of  system  reliability 
it  is  always  assumed  that  the  components  are  tested  separately;  for  a  survey  and 
references  see  Mann,  Schafer,  and  Singpurwalla  (1974).  The  idea  of  basing  the 
estimate  of  the  system  life  distribution  on  estimates  of  the  component  life  dis¬ 
tributions  with  the  resulting  censoring  considerations  is  new  in  reliability  theory. 
This  approach  extends  and  gives  a  novel  application  of  censoring  methodology. 

The  competing  risks  model  corresponds  to  a  series  system.  Aalen  (1976) 

showed  that  for  this  model,  the  vector  of  Kaplar.-Meier  estimates  (F,,...,F  ), 

1  m 

when  normalized,  converges  to  a  multidimensional  Gaussian  process,  whose  com¬ 
ponents  are  independent.  This  result  corresponds  to  our  Theorem  1  for  the  case 
of  a  series  structure. 

This  paper  is  organized  as  follows.  Section  2  gives  some  definitions,  pre¬ 
liminary  results  including  the  strong  consistency  of  F ,  and  results  without  proof 
concerning  the  Kaplan-Meier  estimator  to  be  used  subsequently.  Section  3  is  tech¬ 
nical,  and  applies  martingale  theory  to  obtain  the  results  of  Theorems  1  and  2. 

It  contains  all  the  terminology  and  facts  concerning  martingales  that  are  needed 
to  prove  Theorems  1  and  2.  Section  4  gives  an  application  of  the  results  of  Sec¬ 
tion  3  to  system  design  methods.  The  Appendix  gives  a  proof  of  a  result  used  in 
Sections  2  and  3. 

2.  PRELIMINARIES  AND  DEFINITIONS. 


Corresp  nding  to  a  generic  system,  we  define  generic  random  variables  X  -  ,  2., 

:>  .  ,  and  T,  such  that  the  random  vector  fX  .  ,Z  . ,  <5 .  ,T)  has  the  same  distribution  as 
J  111 


(X; 


i-j  ,Zi  j  *  6i  j  ’Ti )  for  i  =  1>  2,  n>  and  j  =  l>  2,  ...,  m. 


6 


In  Section  1  it  is  noted  that  the  strong  consistency  and  weak  convergence 
results  for  the  Kaplan-Meier  estimator  are  valid  under  the  assumption  that  the 
lifelengths  and  the  censoring  random  variables  are  independent.  In  our  model, 
for  each  j,  X^  is  censored  by  T,  and  for  a  coherent  structure  these  two  random 
variables  are  dependent.  However,  it  is  possible  to  redefine  the  censoring  vari¬ 
ables  to  circumvent  this  difficulty.  This  is  best  explained  in  terms  of  a  simple 
example.  Consider  the  structure  shown  diagrammati ca 1 lv  in  Figure  2.1. 


In  the  example  T  =  X^a  (X,  v  X.)  ,  where  x  a  v  =  min(x,v)  and  x  v  y  =  max(x,v)  .  Consider 
now  component  1.  Clearly  X^  is  censored  by  =  X.,vX,t  which  is  independent  of  X^. 
Similarly,  \?  is  censored  by  Y?=Xj,  and  X,  by  Y  =  X  . 

In  the  Appendix  it  is  shown  that  in  general,  for  each  i  =  1,  ...,  m  there  is 


a  nonnegative  random  variable  Y.  such 


that 


(Z.,6.)  =  (X.  a  Y.  ,  I(X.  <  Y.1), 
y  j  l  j  i  i 


(2.1) 


X.  and  Y.  are  independent.  (2.21 

be  refer  to  Y.  as  the  -jenscring  variable  af  7...  Statements  (2.11  and  (2.21  imply 
that  the  censoring  of  a  component  lifelength  is  described  by  the  random  censorship 
model  (Gilbert,  1962).  Roughly  speaking,  Y  is  the  lifelength  of  the  system  if  X. 
is  replaced  by 

In  order  to  describe  the  distribution  of  Y.  we  introduce  some  notation.  For 

l 

>'  =  f >'  j . ym)  e  f n  ,  1  ] m ,  a  <  [ 0 , 1  ]  ,  and  i  =  1  ,  ....  m ,  let 


Let  F(t)  =  CP, L t ) . F  (t))  and  recall  that  1!.  is  the  distribution  of  I..  In  the 

—  1  m  )  l 

Appendix  it  is  shown  that 

P  (Y  -  >  t)  =  h  ( 1 .  ,  F(t))  ,  (2.4) 

l  f  l  — 

where  h  is  the  reliability  function  (see  (1.4)).  Thus, 

P 

H.(t)  =F.(t)hi(li,  F(t)).  (2.51 

he  now  review  some  terminology  from  reliability  theory  (see  e.g.  Barlow  and 
Proschan,  1981)  to  be  used  in  the  proof  of  consistency  of  F  and  in  later  sections. 
For  a  coherent  system  of  m  components,  the  states  of  the  components  correspond  to 
a  vector  U=  (U^ , . . . ,Um) ,  where  IK  =  I (component  j  is  in  a  functioning  state).  The 
sfructuv s  function  is  defined  by  4>(U)  =  I  (System  functions  when  U  describes  the 
states  of  the  components)  for  U  e  A^,  where  A^={0,l}m.  It  is  well-known  (and  easy 
to  see)  that  for  p =  (p  , . . . ,pm)  e  [0, l]m, 

m  Ii.  1-U. 

h ,  (p)  -  l  0(U)TTp  (1  -  P:)  (2.6) 

Ue A  j  =  l  J  J 

m 

where  0^  =  1  by  definition. 

The  Kaplan-Meier  estimates  F_.  given  by  (1.3)  will  be  denoted  F^  when  we  want 
to  emphasize  the  dependence  on  n;  similarly  for  the  estimate  F  of  system  life  dis¬ 
tribution.  .Also,  Fn  will  denote  the  vector  (F™, . . . , Fn)  . 

—  1  m 

The  following  propositions  specialize  the  properties  of  the  Kaplan-Meier 

estimator  investigated  in  the  literature  to  the  estimators  Fn . 

J 


PROPOSITION  2.1.  (Foldes,  Rejto,  and  Winter.  1980) .  If  T  >  0  is  such  that  II. (T)  <  1, 


then 


ft 


P  :  sup  !  F.  (t)  -  F. ft)  |  =  o  | 
j  0*tt<T  '  -  J  [ 


in  n)  ? 


=  1  . 


Thus,  F .  is  a  strongly  uniformly  consistent  estimator  of  F.  on  the  interval  fO,T|. 


Under  the  assumption  of  continuitv  of  the  distributions  of  the  component 


lit'elengths  and  the  censoring  random  variables,  the  rate  of  convergence  is 
improved . 


PROPOSITION  2.2.  I Foldes  and  Rejto,  1981).  If  F  and  h(l.,Fl  are  continuous  and 
T  >  0  is  such  that  H .  (TO  <  1,  then 

l 

f  r  i 

n  i  1  ;.n  r  „  I  /’ii  Vn  n  I  !  . 

Pi  sup  |F.(t)-F  (t)  =0,  / - - - -i>=l. 

u  0<t<T  '  J  •  i  jj 

We  note  that  if  F  ^ ,  ....  F^  are  continuous,  then  h(l.,F_)  is  continuous:  see  Lemma 
2.1  below. 


PROPOSITION  2.5.  (Breslow  and  Crowley,  1974).  If  F.  and  h(lj,F_)  are  continuous  and 

\*  n 

T  >  0  is  such  that  H.(T)  <  1,  then  n2(F^  -  F^)  converges  weakly  to  a  zero  mean  Gaus¬ 
sian  process  W  on  P[0,T]  whose  covariance  function  is  given  by 


Cov(W(t  ),W(t,))  -  F  (t  )F  (t0)  / 

1  -  J  1  -  n 


dF.fu) 


J  ’  0  H.(u)F.Cu) 

The  following  lemma  is  needed  in  the  proofs  of  strong  uniform  consistency  and 
weak  convergence  of  F. 


LLMMA  2.1.  For  any  structure  of  m  independent  components,  the  corresponding  relia¬ 
bility  function  Ip  is  twice  continuously  differentiable  over  [0,ljU,  and  the  first 
and  second  part ; a  1  derivatives  are  bounded  in  absolute  value  by  1  uniformly  over 

I'M!"'. 


P roo f :  F or  n =  i  p 


,  p  j  c  f  0 , 1  ] "  and  k=  1,  2,  ...,  m,  we  have  by  (2.61, 


— -  ■  =  h  1,  ,n)  -  h  jo  ,n)  . 

P  ■  k  '  k  • 


i  wi:  navi; 


and  for  ..  x  k , 


! 

:.p-~ '  =  -hj'  1k* 1 ;  ’P^  -M  Wp)’-  -  {h5C°k. 1  ,.p)  -  h  - (nk’",  >r! '  -  i:.1.)' 

in  an  ob\-ious  extension  of  the  notation  (2.3).  Bv  (2.6)  is  continuous  over 
m 

[0,1]  .  This  fact  together  with  (2.-),  (2.8),  and  (2.9)  imply  that  the  first  and 
second  partial  derivatives  are  continuous  on  [0,l]m;  hence,  by  Theorem  6. IS  of 
Apostol  (1964),  h^  is  twice  continuously  differentiable  on  [0,l]m.  Equation  (2.-) 
implies  that  the  first  partials  are  bounded  in  absolute  value  by  1.  Since  each 
of  the  two  quantities  inside  the  braces  on  the  right  side  of  (2.9)  is  between  0 
and  1,  it  follows  that  the  second  partials  are  also  bounded  in  absolute  value  bv 
1.  The  lemma  follows  since  p  is  arbitrary.  n 

We  now  establish  the  strong  uniform  consistency  of  F  and  give  the  rate  of 
convergence . 


I’RUPOf  IT  ION  2.4.  Let  T  >  0  be  such  that  for  i  =  1 ,  .  .  .  ,  m,  min(F  .  (T)  ,  F  .  (T))  >  () . 

Then 

f  i  i  r,. 

(a)  P|  sup  |  F(t)  -  F  (t)  |  =  0|  .-n!  —  i  j  =  1 . 

I  0<t<T  1  ,  n '  , j 

(b)  If  F  ...,  F  are  continuous,  the  rate 

f  r  )  ^  _ _ > 

o  j  (  ,n  ,n^ — j  may  be  replaced  by  0  i  ^  -n  —  — —  |  . 

11 "  1  l  n  j 

■••root  :  e  fix  t  |  0 ,  T  ]  and  consider  h  k  (  F_n  1 1  j  I  -  h  ,  ( F_(  t )  1  .  Since  h.  is  conti  minus  1\ 

differentiable  over  [9,1]  by  Lemma  2.1.  we  can  apply  the  '■lean  Value  Theorem  (see 
for  example  Apostol,  1964,  Theorem  6 . IT) :  there  exists  a  point  X*  lying  on  the 
lire  segment  joining  F  ( t 1  and  Fit-  such  that 


in 


where  Vh  is  the  gradient  of  h, .  In  view  of  Lemma  2.1,  Proposition  2.1  proves 
<? 

Part  (a),  and  Proposition  2.2  proves  Part  (b) .  □ 

3.  WEAK  CONVERGENCE  RESULTS. 

In  Section  3.1  we  outline  the  proofs  of  Theorems  1  and  2  and  indicate  where 
the  theory  of  stochastic  integration  and  counting  processes  is  needed.  Section  3 
reviews  the  elements  of  this  theory  that  are  needed  in  this  paper.  Section  3.3 
uses  the  results  stated  in  Section  3.2  to  give  rigorous  proofs  of  Theorems  1  and 
Throughout  Section  3  we  adopt  the  convention  that  jj-  =  0 .  The  index  n  used  in 
defining  a  process  is  suppressed  whenever  possible. 


3.1.  Sketch  of  the  Proofs  of  Theorems  1  and  2. 


To  prove  Theorem  1,  we  show  that  for  any  T>0  satisfying  max  F.(T)  <  1, 

1<  j  <m  -1 


that 


F  -  F 
1  1 


F  -  F 
m  m 


("!•  ••••  V- 


(3.1) 


where  W* ,  ...,  W*  are  independent  mean  zero  Gaussian  processes  with  covariance 
given  by 

tl  dF.(u) 

Cov(W*(t.) ,  b’*(t9))  =  /  - 1 -  forO<t  <t,<T.  (5.2) 

J  J  “  0  H.  (u)F. (u) 

3  J 

(Now  and  henceforth,  the  symbol  d  signifies  weak  convergences  in  Dm[0,T].) 

Theorem  1  is  an  easy  consequence  of  (3.1)  and  Theorem  5.1  of  Billingsley  (1968). 

We  prove  (3.1)  by  a  general  method  introduced  by  Aalen  (1978)  and  later 
refined  by  Gill  (1980).  We  define  the  stopped  process  Ffn  on  [0,*),  j  =  1 ,  2,  ... 
m,  by 


F j  (t) 


F.  (t  a 


1  (n)  i 


and  Ft  ft)  =  1  -  Ftft) , 


(3.3) 


and  use  the  following  decomposition: 


(3.9) 


Nij(t)  =  r(Zij  -  t>  = 


»"(t)  »  I 


i=l  V 


(3.10) 


v..(t)=i(z..,t). 


v"(t)  =  [v..(t). 


1  '  '  n 


(3.11) 

(3.12) 


t 

Pi.fsr 

A.  . 

(t)  =  j 

_ 

dF.  (s)  . 

13 

0  ! 

_Fj(s)  J 

3 

t 

rv"(sr 

f  B  1 

=  / 

_j _ 

dF.(s) 

=  l  A.  At) 

0 

F.(s) 

3 

l  i=i  U  J 

(3.13) 


(3.14) 


M.  . (t)  =  N .  . (t)  -  A.  . (t)  . 

ij  13  13 


f  n  ) 

M?(t)  =N"(t)  -  A?(t)  I  =  I  N.  .  (t) 


i  =  l 


13 


(3.15) 

(3.16) 


The  process  N^(t)  records  the  number  of  uncensored  failures  of  component  j  up  to 


time  t.  The  process  V\  (t)  records  the  number  of  systems  in  which  component  3  is 


at  risk  at  time  t  .  We  note  that  for  each  t,  is  the  completion  of 


a{N^(s);  l<i<n,  l<j<m,  s<t) 


(3.17) 


The  following  proposition  is  fundamental  in  establishing  that 


\  h 


f  F.(t)  -  Ft  (t) 


F*(t) 


Ft  ;  t  e  [0,T] 


>  is  a  martingale. 


PROPOSITION 


3.1.1.  (Gill,  1980).  Suppose  F j ,  ....  F^  are  continuous  and  t  > 0  is 


such  that  max  F.ft)  <  1.  Then  for  each  j  and  for  all  n, 


1<  3  <m 


F.ft)  -  FT  ft) 

-  n’5  r 

Jj(s) 

F  j  ( s  -) 

Ft  ( t) 

-  n  J 

0 

V.  (s) 

F  .  ( s) 

dM-(s) , 


(3.18) 
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where  F.(s-)  is  the  left  limit  of  F.(s). 
J  J 


The  theory  of  counting  processes  is  needed  to  show  that  { (NL  (t)  ,  F^)  ;  te  [0,T]} 
is  a  martingale.  We  apply  the  results  in  Section  3.2  to  show  that 


f 

r 

n'2/ 

0 

V. 

( 

J  j  (s)  F . (s-) 


V.(s)  F  (s) 


dMj(s),  Ft 


;  t  e  [0 ,T] 


>  is  also  a  martingale,  and  to  verify 


the  conditions  of  the  martingale  central  limit  theorem. 

For  fixed  t,  Theorem  1  and  a  standard  application  of  the  delta  method  (see 

1.  ~ 

for  example  Section  6. a. 2  of  Rao,  1973)  yields  the  result  that  n2(F(t)  -F(t))  con¬ 
verges  weakly  to  a  normal  distribution.  We  generalize  this  argument  for  the  pro¬ 
cess  to  prove  Theorem  2. 


3.2.  Review  of  the  Theory  of  Counting  Processes  and  Martingales. 

References  for  the  material  below  are  Chapter  2  of  Gill  (1980)  and  Chapters 
1  -  4  of  Chung  and  Williams  (1983).  A  very  accessible  review  is  Andersen  and 
Borgan  (1984). 

For  a  complete  probability  space  a  family  of  sub-sigma-fields 

{G  ;  t>0}  of  G  is  called  a  standard  filtration  if  for  each  t, 

(i)  (d'.G  , P')  is  a  complete  probability  space, 

(ii)  Gg  c  G  for  s  <  t , 


A  process  Y =  { Y ( t ) ;  t;fl;  is  said  to  be  adapted  to  the  standard  filtration 
•rG^;  t  >  if  for  each  t,  Y(t)  is  -measurable.  The  word  adapted  will  refer  to 
the  filtration  (G  ;  ten). 

A  special  class  of  adapted  processes  is  the  class  of  predictable  processes. 
Rougly  speaking,  Y  is  predictable  if  its  value  at  t  is  determined  by  its  values 
at  times  up  to  but  not  including  f.  \  formal  definition  is  given  on  page  28  ot 
Chung  and  Williams  (19831.  We  will  use  the  fact  that  an  adapted  process  that  has 
a.s.  left  continuous  paths  i>  predict  ible. 


A  martingale  V  is  defined  to  be  square  integrable  if  sup  E  Y“(t)  <  03  • 

t 

Henceforth,  we  assume  that  all  martingales  in  this  section  are  with  respect  to 
(G  ;  t  b 0}  and  are  square  integrable. 

It  is  well  known  that  for  a  martingale  Y,  there  exists  a  unique  predictable 

process  with  nondecreasing  paths,  called  the  quadratic  variation  process  of  Y  and 

•> 

denoted  <Y>,  such  that  -  <Y>  is  a  local  martingale  (see  page  19  of  Chung  and 
Williams,  19S5,  for  a  definition  of  a  local  martingale).  It  follows  from  Propo¬ 
sition  1.8  of  Chung  and  Williams  (1983)  that  uniformly  bounded  local  martingales 
are  martingales.  Since  all  local  martingales  encountered  in  Section  3.3  are  uni¬ 
formly  bounded,  the  reader  may  substitute  the  word  martingale  for  local  martingale 
throughout  this  section  without  affecting  the  material  in  Section  3.3. 

If  Y  and  Y.,  are  martingales  then  their  covariation  process  <Y  ,Y?>  is  defined 
by  <  Y  j ,  Y  ,  >  =  h  ( <  Y  j  +  Y .,  >  -  <  Y  j  -  Y .,  > )  .  Kunita  and  W'atanabe  (1967)  showed  that  sY^.Y^s 
is  the  unique  predictable  process  with  paths  of  bounded  variation,  such  that  the 
process  Y^Y,-<Y^,Y^>  is  a  local  martingale.  It  is  easily  seen  that  for  every 
martingale  Y,  we  have  <Y>  =  <Y,Y>. 

For  a  process  of  bounded  variation  Y,  define  ||y|(  to  be  the  process  such  that 
' ’  Y ; !  ( t )  is  the  variation  of  the  paths  of  Y  on  [0,t].  The  following  proposition  is 
used  to  show  that  the  integral  in  (3.18)  is  a  martingale. 


PROPOSITION  3.2.1.  (Doleans-Dade  and  Meyer,  1 9 “'0 )  .  Let  G '  and  G'  be  predictable 
processes  and  Y  and  Y.,  be  square  integrable  martingales  of  bounded  variation  such 


that 


J  j  G'(s)  jd|!  Yk||  (s) 


<  <=»,  k  =  1 ,  2 . 


Then  the  processes  and  defined  on  [0,t]  whose  paths  are  defined  by  the 
Lebesgue-Stielt jes  integral 


X,  (s)  =  /G'(u)dY.  (u)  s-t, 


is  a  local  martingale  on  f 0 , t ] .  Furthermore, 
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s 

<X  ,X9>(s)  =  /  Gl'(u)G'(u)d<V  ,  Y .,  >  (  u  )  , 

0  “ 

for  s  <  t . 

A  vector  of  adapted  processes  (Y^,...,Y  )  is  called  3  counting  process  if  the 
following  hold  a.s. 


(i) 

Y.  (0)  =0,  i  =  1,  2 

,  .  ■  •  ,  m . 

(3.19) 

(ii) 

The 

paths  of  each 

process  Y.  are  nondecreasing, 

right  continuous 

and 

have  jumps  of 

sice  +1  only. 

(3.20) 

(iii) 

No 

two  processes 

jump  at  the  same  time. 

(3.21) 

Theorem  1.9  of  Meyer  ( 1 9 T’fs )  implies  that  for  each  process  Y.,  there  is  a  unique 
predictable  process  B.  with  right  continuous  and  nondecreasing  paths  originating 
at  0,  such  that  Y.  -  Bj  is  a  local  martingale.  The  process  B.  is  called  the  com¬ 
pensator  of  Yj  . 

The  following  proposition  is  adapted  from  a  theorem  of  Murali-Rao  (1969)  and 
is  useful  for  identifying  the  compensator. 

PROPOSITION  5,2.2.  (Gill,  1980).  Let  Y  be  a  univariate  counting  process  and  let 
t  e  (0,=°)  satisfy  E  Y(t)  <°°.  Define 

t,  =  22‘kt,  k  =  1 ,  2.  ....  2  =  0,  ],  ....  2k 

anJ  2k-I 

V  I  E(T(tk.  ...I'’  Yltk,t>  !Gtt  >•  k*>'  2 . 

Then  there  exists  a  subsequence  of  integers  fr^l,  as  k-*-T’,  and  a  unique  ran¬ 

dom  variable  l),  such  that  for  all  bounded  random  variables  X, 

F.  (XU  )  -*■  F.  ( XU ) 

rk 

as  k  -  .  The  compensator  B  of  Y  satisfies 

B(t)  =  II  a.s.  . 

We  note  that  this  result  holds  for  each  fixed  value  of  t,  and  that  special 
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care  needs  to  he  taken  when  dealing  with  a  continuum  of  t's. 

The  following  proposition  is  a  special  case  of  Theorem  2.3.1  of  Gill  (1980). 

PROPOSITION  3.2.5.  (Gill,  1980).  Suppose  (Y  , . . . .Y^)  is  a  counting  process  with 
compensators  (B^ , . . . , 8^) .  Define  the  process  Z  .  by  Z .  = -  B . ,  i  =  1 ,  2 ,  . . . ,  m. 
If  the  processes  (B  ,...,Bm)  have  a.s.  continuous  paths  then  the  following  hold. 

(i)  Z.  is  a  local  square  integrable  martingale,  j  =  l,  2,  ...,  m. 


tii)  <Z.  >  3  >  = 

J 1  -  2 


f  V  ifjl  =  j2 


0  if  ji  * j2 


We  use  a  multivariate  extension  of  a  martingale  central  limit  theorem  due  to 
Rebolledo  (1980)  in  Section  3.3.  The  conditions  of  the  proposition  below  are 
stronger  than  those  used  by  Rebolledo.  We  use  the  stronger  conditions  because 
they  are  easier  to  understand  and  do  not  take  much  effort  to  verify  in  Section  3.3. 
Let  f  ,  ....  fm  be  positive  functions  on  [0,®). 

PROPOSITION  5.2.4.  (Gill,  1980).  Suppose  that  the  sequence  of  vector  processes 

(Zj, . . . ,  Z^)  n  =  1,  2,  ....  satisfy  the  following  conditions .  For  every  t  >  0, 

1  s  j  ,  j,5m,  s  e  [0, t] ,  and  every  n , 


(i)  Z.  is  a  square  integrable  martingale, 


....  ,n  _n  .  .  p  fj(5)  if  jl*G 
(li)  <Z.  ,  u.  >  (s)  U.  < 


J,  y 


(3.22) 


(3.23) 


1  f  -’l^; 


there  exists  processes  z”£,  Z*?£  such  that 

(iii)  Z*?£  and  Z_j£  are  square  integrable  martingales, 

, .  ,  „n  =ne  _ne 
I  1  V)  Z  •  =  ~  .  +  L  ■  , 

J  J  -J  ’ 

(v)  has  no  jumps  larger  than  e, 

(vi)  ZjE  has  a.s.  paths  that  are  of  bounded  variation, 


(3.24) 


(3.25) 


(3.26) 


(3.2-) 
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(vii)  the  processes  Zj£  and  Z?e  do  not  jump  at  the  same  time, 


,  •  •  ■  =ne  -ne  ,  .  p  _ 

(vm)  >(s)  i  0,  as  n 


(3.28) 


(3.29) 


Then, 


n  n  ti  °°  00 

(z;,...,z^)  *  (Zr...,zm) 


ou  00 

where  Z ^ ,  ...,  Zm  are  independent  Gaussian  processes  with  mean  zero,  and  covariance 


structure  given  by 


Cov(Zj  (sp  ,  Z.  (s2))  =  f  .  (Sj)  0<Sj<s2<t. 


3.3.  Weak  Convergence  Results. 

All  families  of  o-fields  defined  in  this  section  are  standard  filtrations, 
This  fact  is  a  consequence  of  Theorem  A. 2.1  of  Gill  (1980). 


LEMMA  3.3.1.  Suppose  that  F ^ ,  ....  F^  are  continuous  and  suppose  T  is  such  that 


Fj (T)  <  1  f or  j  =  1 ,  2,  . . . ,  m .  Then  for  each  n  «  1 ,  2 ,  ...»  we  have 

fl+MT) 


(i)  The  variation  of  the  paths  of  M11  is  bounded  by  m 


F  •  (T) 


on  [0  ,T]  .(3.30) 


(ii)  {(Mj(t),  F^) ;  te  [0,T]}  is  a  martingale  for  each  j,  where  F^  is  defined 


by  (3.17) 


(iii)  <M"  (t),  Mj  (t): 
J 1  J2 


,  J  A"j(t)  = 


if  h*  j2 


(3.31) 


(3.32) 


Proof :  The  proof  of  (i)  is  immediate. 

To  prove  (ii)  ,  we  show  that  for  i  =  l,  2,  ....  n,  NT  is  a  martingale  on  [0,T]. 
Part  (ii)  then  follows  since  a  sum  of  martingales  is  a  martingale.  If  we  assume 
that  Fj  is  absolutely  continuous  with  a  continuous  derivative  and  that  the  compen¬ 
sator  of  Nhj  is  absolutely  continuous  with  a  left  continuous  derivative  whose  right 
limit  exists,  a  simple  proof  that  NT  .  is  a  martingale  on  [0,T]  can  be  given  as 


fol lows . 


It  is  not  difficult  to  verify  for  every  t?  [0,T]  that 


.vr.’-v:  j-  .* 


aJK* .  “A. .  L  . . 


IS 


1  fi(t) 

1  im r- P (N .  .  (t  +  h)  -  N.  (t)  >  l|  F  )  =  V  (t)  -* -  , 

h-'O  -  -  -  F.(t) 

where  f.  denotes  the  density  of  F . .  This  is  enough  to  show  that  A. .  is  the  com- 
1  i  '  il 

pensator  of  Nf  .  .  (See  Lemma  3.3  of  Aalen,  1978.')  It  now  follows  that  is  a 
martingale. 

We  now  consider  the  general  case.  For  measurability  reasons  that  are  indi¬ 
cated  later  on  in  the  proof,  we  first  show  that  is  ?  martingale  with  respect 
to  the  filtration  {G  ;  te  [0  ,T] } ,  where  for  each  te  [0,T],  G^  is  the  completion  of 
a{I(X^j<s);  1  <  i  <  n,  i  <  j  <  m,  s  <  t)  .  We  then  use  this  fact  in  a  simple  argument 
to  show  that  Nb  .  is  a  martingale  with  respect  to  (F^;  te  [0,T]}. 

To  prove  that  is  a  martingale  with  respect  to  (G^;  te  [0,T]},  we  show 
that  A„  is  the  compensator  of  N^.  with  respect  to  (G^;  te  [0,T]}.  Let  denote 
the  compensator  of  N.^  with  respect  to  (G^;  te  [0,T]}.  The  key  step  in  proving 
that  A..  =U..  a.s.  involves  the  use  of  Proposition  3.2.2  to  show  for  anv  fixed 

ij  ij  K 

t  e  [0,T]  that 


PCA.jCt)  =u..(t))  =  1. 


(3.33) 


Since  A.,  has  a.s.  continuous  paths,  it  follows  that  A..=U..  a.s.  on  [0,T1.  To 
ij  r  >  ij  ij  l  >  j 

prove  (3.53)  we  use  the  following  definitions.  For  te  [0,T]  define 

tj.  t=t2"kt,  k  =  1 ,  2,  ....  1=1,  2.  ....  2k, 


v  ;„E‘Kij<tk.t.1)-Nu(tk,AiG,  >  k=>,  2 . 


A  key  step  in  proving  (3.33)  is  to  show  that 

U,  5  A. .(t),  as  k 

k  ij 


(3.34) 


Assume  that  this  has  been  done.  Proposition  3.2.2  then  implies  that  there  exists 


a  subsequence  of  integers  (r^),  r^  -*■  ®  as  k  -*■  ■*>  and  a  random  variable  V  satisfying 


E  (XU  )  ->■  E  (XV)  ,  as  k  -*  =° 
rk 


(3.35) 


for  all  bounded  random  variables  X,  and 


pCv  =  Uj . (t) )  =  l. 


(3.56) 
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It  follows  from  (5.34)  and  (3.35)  by  standard  probability  arguments  given 

below  that  A. . (t)  =V  a.s..  Statement  (3.33)  now  follows  from  (3.36). 
ij 

For  the  sake  of  completeness  we  show  that  (3.34)  and  (3.35)  imply  that  (t) 
V  a.s..  It  follows  from  (3.34)  and  (3.35)  that  there  exists  a  subsequence  of 


integers  {A„},  as  l -*■<*>  such  that 


U  A  (t)  a.s.  , 


(3.37) 


E  XU.  •+  E  XV, 

Az 


(3.38) 


as  Z  ■* °°  for  all  bounded  random  variables  X.  Let  e  >  0  be  arbitrary.  Define  the 


sets  CjCC,c  as  follows: 


C  =  n  {|U,  -A. . (t) |  <  e),  a  =  1,  2 . 

a  l=ct  AZ  1J 


(3.39) 


Statement  (3.37)  implies  that 


lim  P(C  )  =  1. 
a 

a-*® 


(3.40) 


Let  D  and  D  denote  the  sets  {V>A^.(t)  +  e}  and  {V  <  A^.  (t)  -  e }  ,  respectively. 
For  each  a,  (3.58)  implies  that 


lim  E (I  (D+  n  Ca)UA  )  =  E  I(D  +  n  C  )V. 
a  A  l 


(3.41) 


Since  the  definition  of  implies  that 


lim  E(I(D  nC)UA  )<E(I(D  n  C)  (A .  .  (t)  +  e) )  , 
L -KB  l 


it  follows  from  (3.41)  that 


E ( I  (D  nCJV)  <E(I(D  n  C^)  ^  •  (t)  ♦  c) )  . 


(3.42) 


We  conclude  from  the  definition  of  D  and  (3.42)  that  P(P  n C  )  =0  for  each  a. 
Statement  (3.40)  implies  that  P(D  +  )  =0.  Similar  arguments  show  that  P(D  )  =0. 


Since  c  was  arbitrary,  we  conclude  that  Y  =  A..ft)  a.s. 


We  now  prove  (3.34).  For  each  k=  1,  2,  ....  i=  1,  2,  ....  detine  the  set 

Bl  +1  =  >  tl  -  +  1  ]  .  ;  >  t.  . 

K  11  K  »  *-  K  »  *  1  11  k  ,  ■ 


and  the  randor.  variable 


>'  (xijltk.;*i)-su(tk.d) 


We  note  that 


’k ,  1 


i  1  if  t,  .  <  V.  .  <  X.  .  <  t 
=  J  lj  li  k,.+l 

10  otherwise 


(3.43) 


(where  if  .  is  defined  by  (2.1)  and  (2.2))  and  that 


"  ,  is  G  -measurable  for  all  k  and  l . 

k ,  <-  t,  , 

k ,  c 


(3.44) 


It  is  ease  to  see  that  ,  is  not  F  -measurable  and  this  is  the  reason  we  use 

K  , t,  , 

k  ,  >. 


the  filtration  {G  ;  te  [0,T]}.  To  prove  (3.34)  we  first 


show  that 


oL  "k>  i+i 


0 ,  as  k  -*■ 


(3.45) 


Assuming  the  validity  of  (3.45),  Lemma  2.5  of  Rootzen  (1983)  together  with  (3.45: 
and  3.44  immediately  imply  that 


Fi  (1(5 


k  t 


k ,  j 


-xii(tk,dlGt.  .'j Eo'  « 

=1  •  '  k , Z 


k^<*>.  (5.40) 


We  then  show  that 


A .  .  !  t  i  -  ;  F.  (  I  (  B,  )  G  )  •  -  0  a  .  s  . 

u  ;  =n  k  rk , 


(3.4-) 


Statements  (3.46)  and  (3.4"')  now  imply  (3.34). 


We  now  prove  (3.45).  For  each  k 


-  21  - 


2k-l  tk,i+l 


i=  0  t.  , 

k , 


/  (F  (t  )  -  F  (u))dh  (1  F(ul) 

J  J  K  ,  1  J  -5  1  ~ 


I  [F.(t.  .  ,)-F.(t,  J][h  fl.t  F  (t,  J)-h  (1  ,  F(t  ))] 
„L  n  l  k,i  +  l  i  k,^  5  i  —  k,t  t>  i  k,«.+  l 

i= 0 


sup  (F  .  (s  +  -T-)  -  F  .  (s) )  -*•  0  as  k 
l  1 


S  6  [  0  ,  t  -  -£] 


since  F.  is  uniformly  continuous  on  [0,t]. 

We  now  prove  (5.47).  Since  for  each  k  and  2.  I (Z -  -  >  t.  )  is  (  -measurable, 

aj  K.i  ^k , Z 

we  have 


It  is  easv  to  verify  that 


tk ,  £  +  1  dF.(u) 


k ,l  Fi(tk,i') 


Thus , 


E(I<\  >l\  ’’"kj-k,.1  I 

k ,  Z  t. 


k,2.+  l  dF  .  (u) 
r  J 


h,l  Fj(tk,t) 


Define  the  random  variables  Y,  by 

k 


2-1  k, 2+1  dF.(u) 

V  E  Hzu  ^ — 

K  1=0  J  ’  t.  .  F.(ui 

k  ,  i  1 


We  prove  (3.47)  by  showing  that 


l  E(I  (Bf'  )  |G  )  -  Y.-O  a .  s .  as  k 

2=0  k  \,l  K 


Y,  -+A.  .  (t)  a.s.  as  k-*°°. 
k  ij 


We  proceed  to  prove  these  assertions. 


- 1  k,  .1  +  1 


n  1 1 '  b. 


a!  .it; 

i 


1  1 

<  sup  - T—  -  - - 

0<s<t  -  ~  !  F  .  I's  +  F.(V 

k  —  i  _,k  j 


0  as  k  -* 


the  unitor:: 


continuity  of  F.  on  [ 0  ,  t ] .  Next, 

j  _  _J_ 

1  '  *  ,k  d;  .  (ui 

t  ,  .  , .  1  j  1  I 

:  A .  .  it )  -  i .  :  <  sup  ,  j  — —  1 

11  k  ‘  1  i  J  r  f  \  \ 

0<s<t - r-  I  S  F  .  (u) 

->K  J 


=  — - —  sup  (F  .  ( s  +  -~y)  -  F  .  ( s ) )  -  0  as  k  -  ^ . 

F  •  ( t )  0<s<t  — r-  J  2  J 

J  -»* 

rhis  proves  to.-!"),  consequently  (3.54),  and  hence  (5.53). 

he  now  prove  that  is  a  martingale  with  respect  to  the  filtration 
■  F  ;  t  [P,T]'.  It  is  easy  to  \rerify  that 


F  = G  for  each  t  e  [0 ,T] , 


( 3 . 4  S 


N.  .  is  adapted  to  ;F^;  te  [0,T]V 


(3.49 


anc  tnai 


V..  is  adapted  to  rF  ;  te  fn.T]''. 


(3.50 


It  t  ;  1  1  'ws  rr  -  '3.49)  and  i  3 . 50  .<  that 


M.  .  is  adapted  to  ;F  ;  te  Ul.Tp 
i  ’  ■  t 


13.51 


It  now  follows  from  (3.48)  and  (5.51)  that  for  each  i  and  all  set  v  T  that 

F  (M .  .  ( t )  '  F  ,  =  F.  f  E  f  M .  .  C  t )  !  G  ) 1  F  >  =  M .  .  i  s )  . 

11  S  1  1  S  S  1  J 

Tli  us ,  M.  .  i  s  *  ".art  ingale  with  respect  to  <T  •  te  ( 0  ,T)  ,  and  Part  (iii)  has  been 


pro’-  cu . 


c- n-L-au-rcc  of  (ii)  is  that  A1?  is  the  conpcnsator  of  Nn  for  each 


toil  ,.s  from  [’reposition  3.2.3. 


An  alternative  way  to  prove  that  M  is  a  martingale  is  to  use  Lemma  3.1.1 

ot  Gill  l  19 SOI  to  show  that  A.  .  is  the  compensator  of  N.  ..  While  shorter,  this 

1  j  lj 

approach  is  not  as  straiehtforward  as  the  method  used  above. 

We  use  the  filtration  (F  ;  te  f(l.T|}  because  it  has  a  natural  interpretation: 
represents  all  the  information  available  at  time  t.  It  is  also  the  filtration 
most  widely  used  in  the  literature.  However  it  is  clear  that  Lemma  3.3.1  is  valid 
for  any  filtration  satisfying  (3.48),  (5.49),  and  (5.50),  and  furthermore  the 
processes  defined  by  (3.S)-(5.16)  are  adapted  to  this  filtration.  For  the  remainder 
of  the  proofs  of  this  section,  the  filtration  enters  only  via  Lemma  5.5.1  and  the 
measurability  of  the  processes  defined  by  (3 . 8) - (3 . 16) .  Thus,  the  results  of  this 
section  can  be  proved  using  any  filtration  satisfying  (3.48),  (3.49),  and  (5.50). 
However,  for  the  reasons  stated  above  we  shall  continue  to  use  { F^_ ;  te  [0,T]}. 


Proposition  5.2.1  together  with  (5.55)  and  (5.18)  imply  that  for  each  n  and  j. 


i  i  . 


I  p* 


is  a  local  martingale  on  [0,1],  Since  these  local  martingales  are  uni- 


J 


formly  bounded  for  each  n,  they  are  square  int eg  table  mart  ingales  on  [0,T],  The 
proof  of  [5.6)  follows  from  Proposition  5.5.4,  whose  conditions  we  now  verify.  We 

,  f  -  r+  , 

1  f  r  :  ,  l 


have  shown  that  (5.22)  holds  for  n  ^  ,  and  we  new  check  (5.25).  Proposi 


tion  5.2.1  implies  that  for  l<i  ,  j^sm, 


F .  -  Ft 


<  n 


•'  1  J 1 


Ft 


.  n 


F  .  -  F* 

}->  J! 


I  2 1  )  l  J  -i  ) 


=  n  j 


r  I  1  ‘  1  i 


(5.54) 


0  |  V  (s ) F  .  (s) 

-  -1  J 1 


J,  (s) F  (s-)' 


J 


V.  (s)F  (s) 
-  -  2  -  2 


d  <  M.  ,  M.  >  (s) 


J1 


Thus,  it  follows  from  (3.32)  that  for  j  *  j_, 

1  u 


f  - 

•  a 

F  .  -  Ft 

F  .  -  Ft 

>1  J1 

jo  j  a 

L.  \ 

Ft 

F  *  1 

r  .  1 

1  ^ 

b  ; 

(t)  =0. 


(3.55) 


Statements  (3.54)  and  (3.32)  yield  for  each  te  [0,T]  that 


,  i  F  .  -  F*  ! 

n  2  i  I  ,  n  2 


Fi  > 


F  .  -  F*  1 
J  ^  I 


Ft 

1 


ft)  =  n  j 


■  I  fs)F  j  (s-) 


0  Mi .  (S) F  .  (s) 


d,\  (s) 


5 . 5b  ) 


n  / 
0 


J . (si rf .  is-)  v 


_v. (s) ( F  .  ( s ) ) 


dF  (s)  . 


The  uniform  consistency  of  the  Kaplan-Meier  estimator  (Proposition  2.2)  and  an 
application  of  the  G1 i venko-Cantel 1 i  Theorem  (modified  in  a  minor  way)  together 


with  the  fact  that  H .  <  1 1  >0  now  give  that 
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t[~J .  (s)F"(s-)l  t  dF  (s) 

n  /  — - 1 -  dF .  (s)  -*■  f  - 2 - 

oLv.(s)Fj(s)  J  J  0  H.  (s)F.(s) 


a.s.  as  n 


Thus,  condition  (5.25)  follows  from  (5.55)  and  (5.5"). 
For  each  e>0  define  the  following  processes. 


U  t  rj.(s)F.(s-n  f  n2J,(s)F.(s-) 


n  2Zn£  (t)  =  n  2  /  -1 


o  Lv-j  Cs)Fi  Cs) 


— * - *! - >  e  d.M  .  f  s '}  . 

Vj(s)F.(s)  j  J 


nVE(t)  =  n  2  / 


t  rj.(s)F,(s-)l  (  n  2 J . ( s ) F . ( s - ) 


Vj(s)Fj(s)  _ j  (  V  .  (s)Fj (s 


<  e  dM  .  (s)  . 
J 


It  is  clear  that  Z^£  and  are  square  integrable  martingales,  and  that  their 
sum  is  equal  to  the  right  side  of  (5.i8).  Conditions  (5.26),  (5.27),  and  (5.28) 
are  trivially  satisfied. 

We  now  check  (5.29).  Propositions  5.2.1  and  (5.52)  imply  that  for  each 


t  e  [0,T] , 


-ne  =ne  .  .  r  I  J 
<Z  ,  -  >(t)  =  n  J  -J -  73 

J  3  0  V. (s)  F 


tfJ.(s)  F‘(s-)1  f  .  J.(s)  F.(s-) 


(s)  F  (s) 


Vj (s)  F  (s) 


•>  e  dF. (s) 


(5.58) 


Almost  surely,  the  indicator  inside  the  integral  is  0  for  all  large  n,  by  the 
G1 ivenko-Cantel  li  Theorem.  This  proves  (5.29)  and  concludes  the  proof  of  Theorem 


Originally,  we  proved  the  asymptotic  normality  of  the  vector  (F^  -F^,  ..., 

F  -  F  )  using  the  method  of  Breslow  and  Crowley  (1974).  The  proof  was  conceptually 
simpler,  not  requiring  the  introduction  of  various  families  of  o-fields  and  the 
heavy  machinery  of  stochastic  integration  and  martingale  central  limit  theorems. 
However,  we  were  unable  to  obtain  the  covariance  terms  in  the  asymptotic  covari¬ 
ance  matrix. 

Suppose  the  life  distributions  of  the  components  all  have  infinite  supports. 

It  is  straightforward  to  show  that  Theorem  1  is  equivalent  to 
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n  2(F,  -  F  >  •  •  •  ,F  -  F  )  -►  (W  ,  .  .  .  ,W  )  weakly  on  D  [0,“)  , 
1  1  m  m  1  m  -  1  J 


(3.59) 


where  D[0,°°)  has  been  equipped  with  the  standard  metric  for  convergence  on  com- 

pacta  (see  Definition  1,  page  123  of  Pollard,  1984).  It  is  not  hard  to  prove  that 

a  functional  f  defined  on  Dm[0,°°)  is  continuous  with  respect  to  the  above  metric 

if  and  only  if  f  is  continuous  with  respect  to  Dm[0,T],  for  each  T>0.  Thus,  (3.59) 

offers  no  advantage  over  Theorem  1  in  obtaining  via  Theorem  5.1  of  Billingsley 

(1968),  the  asymptotic  distributions  of  functionals  of  F, ,  ...,  F  . 

l  m 


Proof  of  Theorem  2:  The  uniform  bound  for  the  first  two  partial  derivatives  of  h^ 
given  by  Lemma  2.1  together  with  Taylor's  Theorem  imply  that  for  each  tc  [0,T], 


m 

3ha 

F(t)  -  F (t)  -  l 

j  =  l 

—A 

3u  (u  ,. 

.  J  _  ^ 

•*V; 

(F.(t)  -  F.(t)) 

(F1(t),...,Fm(t)) 


(3.60) 


km  m 


l  l  sup  F  (t)  -F  (t) 
jj*  1  j2«l  V0<t<T  J1 


sup  F  (t)  -  F  (t) 
0<t<T  J2  J2 


It  follows  from  Proposition  2.2  that  the  right  side  of  (3.60)  converges  to  0 
a.s..  We  use  the  fact  that  convergence  in  sup  norm  implies  convergence  in  the 
Skorohod  topology  (see  page  111  of  Billingsley,  1968)  to  conclude  that  the  process 


!  m  3h 

F  (t)  -  F  (t)  l  — ^  (  (F(t)-F(t))  -*-0  a.s.  in  D[0,T]  . 

j  =  1  j  (ul’*--’um}_  J  J 

(F1(t),...,Fm(t)) 

Thus,  the  proof  follows  by  showing  that 

!  rn  9h.  )  .  , 

n"2  I  (u  u  )  -  (Fi(t) 

j  =  l(3uj  (U1 . Um  J  J 

(Fj(t),...,Fm(t)) 

which  is  the  consequence  of  Theorem  5.1  of  Billingsley  (1968)  and  Theorem  1.  r 
To  construct  confidence  intervals  for  F(t),  we  define  the  following  functions 


(F.(t)  -  F.(t))  XW, 


and  processes  on  [0,™). 


1 


r*  * 
% 


:• 


b 


ra 


t: 


G  (t)  =  (F  (t))“  / 

J  J  rv 


0  F.(s)H.(s) 


t  dN  (s) 

Gj(t)  -  /  3 - *4 - I  =  "  (F,(t)T  l 


0  H.(s-)H.(s) 
J  J 


(i)j  1 


V  "  L  <t(n-  i  +  l)Cn-  i)  J  ' 

mr 


where 


1  n 

H.  (t)  =  -  V  I(Z.  .  >  t)  ; 

J  n  ij 


tyt)  = 


h"(t)  = 


\ 

8u . 

.  J 

(V-.. 

•Um}  = 

(Fj(t), 

...,Fm(t)) 

K 

3u. 

(v--* 

,u  )  = 

(fj(t). 

The  quantity 


G.(t) 

— - -  is  called  Greenwood's  estimator  of  the  variance  of  F . (t)  . 

n  J 


LEMMA  5.3.2.  Suppose  F^,  .  ..,  F^  are  continuous  and  T > 0  is  such  that 


m 


max  F  (T)  <  1.  Then  £  h.G.  is  a  strongly  consistent  estimator  of  [  h.G.. 
1< j  <m  J  j  =  l  J  J  '  j  =  l  J  J 


We  note  that  in  view  of  (1.7),  Theorem  2  and  Lemma  3.3.2  allow  the  formation 
of  asymptotic  confidence  intervals  for  F(t),  te  [0,T] . 


Proof:  Part  (a)  of  the  proposition  in  Section  2  of  Hall  and  Wellner  (1980)  together 


with  Proposition  2.2  imply  that  G.  is  a  strongly  consistent  estimator  of  G_.  .  Lemma 


2.1  implies  that  the  partial  derivatives  of  h  are  continuous.  Thus,  it  follows 


from  Proposition  2.2  that  is  a  strongly  consistent  estimator  of  h...  The  proof 


fol lows  . 


4.  ESTIMATION  OF  THE  RELIABILITY  IMPORTANCE  OF  COMPONENTS. 


The  rel 


eligibility  importance  I  .  (t)  of  component  j  at  time  t  is  defined  by 


I  .  (t)  =^-  ,  ,  '  .  (4.1) 

1  ;U.  hlU, )  I  , 

i  1  m  '  (u, . u  i  = 

i  1  m 

'  ( F.  (tl ,  .  -  .  , F  1 1  i  i 
1  m 

Let  e  be  small  numbers.  Note  that 

1  m 

h^+£l . . frU),;  ]  t,. 

1  =  1  ■  ' 

Thus,  the  reliability  importance  of  components  may  be  used  to  evaluate  the  effect 
of  an  improvement  in  component  reliability  on  system  reliability,  and  can  there¬ 
fore  by  very  useful  in  system  analysis  in  determining  those  components  on  which 
additional  research  can  be  most  profitably  expended.  For  details,  see  pages  26- 
28  of  Barlow  and  Proschan  (1981),  and  the  review  by  Natvig  (1984). 

We  estimate  I  (t)  by  replacing  (F^t) . F^t))  with  (F^t)  ,  .  .  .  ,Fm(t) )  in 

(4.1).  Formally,  define  I _.  by 


L(t>  =  3t-h<V  j  (Uj . 

1  (F1(t) . Fm(t)) 


(4.2) 


PROPOSITION  4.1.  Suppose  Fj,  ....  F^  are  continuous  and  T>0  is  such  that  F .  (T)  <  1, 
i  =  l,  2 ,  ....  m.  Then 

in’ll,  -  I ......  I  -  I  I  -  (Y. . V  ) , 

11mm  1  m 

where  lY, . V  1  is  a  vector  of  mean  zero  Gaussian  processes  whose  covariance 

1  m 

structure  is  given  by 

CovfY.  (t  ),  Y.  (t,)) 

”1  J  ">  *“ 


k-i r,ui  a  '  = 

fVV . 


I  l~h  | 

1  HU  .  3u.  I  (u  ,  .  .  .  ,U  ) 
(  k  !  _  1  m 


f  F  ( t ) . F  It, 

12  m2 


1  dp  (u) 

1  1 1  )  F  ( t ,  I  '  - - - -  ,  for  net  e  t  t  T  and  i  ,  .i  ,  =  I ,  .... 

K  1  k  "  0  f  u)  f  k  ( u )  " 
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As  before,  the  covariance  terms  in  (4.3)  can  be  estimated  consistently,  enabling 
the  construction  of  confidence  intervals  for  I.(t). 


Proof:  For  each  j  define  the  function  g.  bv 


MP1 . ^  *$<“i . V 


(Ujm.mUJ  for  0<pk<l,  k=l,  2,  m.  (4.4) 

=  (Pr---.Pm) 


Assume  that  we  can  show  that 


gj  is  twice  continuously  differentiable  on  [0,l]m  with  first 

and  second  partials  bounded  in  absolute  value  by  1  uniformly  (4.5) 

over  [0,l]m, 


(cf.  Lemma  2.1).  The  proposition  then  follows  by  a  straightforward  multivariate 

” J  '  ■  • • ,  m. 

We  now  prove  (4.5).  It  follows  from  (2.7)  that 


extension  of  the  proof  of  Theorem  2  with  h  replaced  by  ,  j  =  1, 


9gj 

3p7 


d2h. 


and  that 


3‘g, 


9pk3p5 


3p.3p. 
P  k 


33h, 


3p.3p,3p„  ’ 

p  p 


(4.6) 


(4.7) 


for  each  pe  [0,l]m.  It  follows  from  Lemma  2.1  that  g_.  is  continuously  differen¬ 
tiable  and  its  first  partial  derivatives  are  uniformly  bounded  in  absolute  value 
by  1  on  [0,l]m.  Statement  (2.9)  implies  that  for  distinct  indices  j,  k,  f. , 


3\ 


3p  .  3p.  3p  j 

'  .1  rk  tip 


=  VVVVp^VVVV^'VVVVp) 

(4.8) 


fht(0j’1k’V'PJ  '  W  Wp>>  +  {h*(0j»°k'1l»p)  '  WV°*’P)1, 


in  an  obvious  extension  of  the  notation  (2.3).  Statement  (2.8)  implies  that  if 
at  least  two  of  the  indices  j,  k,  i  are  equal  then 


30 


aV 


Vkp* 


=  0. 


(4.9) 


It  follows  in  the  same  manner  as  in  the  proof  of  Lemma  2.1  with  (4.8)  and  (4.9) 
replacing  (2.7)  and  (2.8),  that  h^  has  a  continuous  third  derivative  on  [0,l]m  and 
the  third  partial  derivatives  are  bounded  in  absolute  value  by  1  uniformly  over 


[0,1]  .  Thus,  for  each  j,  has  a  continuous  second  derivative  and  the  second 
partial  derivatives  are  bounded  in  absolute  value  by  1  uniformly  over  [ 0 , 1 ] m - 


□ 


APPENDIX:  RANDOM  CENSORSHIP. 

In  Section  2  we  assumed  the  existence  of  a  censoring  random  variable  Y ^  that 
satisfies  (2.1),  (2.2),  and  (2.4).  Here  we  define  Y.  and  formally  prove  that  it 
satisfies  (2.1),  (2.2),  and  (2.4).  Define  the  binary  function  #.  by 

um)  =  <f|(lj»Uj,...,  um)  ,  u^  =  0,  1,  k  =  1 ,  2,  ...,  m,  (Al) 

where  $  is  the  structure  function.  (See  the  paragraph  preceding  equation  (2.6).) 


The  censoring  random  variable  Y^  is  defined  as  follows: 


Yij=sup{t:  <t>.  (I(Xil  >  t)  ,  ...,  I(X.m>t))  =  1}, 


(A2) 


PROPOSITION  A . 1 .  For  each  j,  Y^ ,  Y  ,  ....  are  i.i.d.  random  variables  satisfying 
(2.1),  (2.2),  and  (2.4). 


Proof:  It  follows  from  (A2)  that  Y. .  is  a  function  of  the  vector 

-  ‘  ii 

( 1  .  ,  I  (X .  ,  >  t ) ,  . . . ,  I (X .  >  1 1 )  .  Thus  it  fol lows  that  Y  ,  Y  ,  . . . ,  are  i.i.d. 

j  i  1  i  m  1  j  -  .1 

and  that  Y..  satisfies  (2.2). 

U 

We  proceed  to  prove  (2.4).  The  structure  function  4>  is  increasing  in  its 
arguments  (see  Definition  2.1,  page  6  of  Barlow  and  Proschan,  1981)  and  hence  a 
fortiori  <Jk  is  increasing  in  its  arguments.  Thus 


P(Y.  .  >  t)  * P(*.  (I(X. j  > t) ,  ....  I (X .m  >  t ) )  =  1 ) . 


(A  31 


It  is  easy  to  see  that  the  right  side  of  (A3)  is  equal  to  h ( 1  .  ,  £(t))  and  so  Y. 


satisfies  (2.4).  To  prove  that  Y_  satisfies  (2.1),  we  consider  two  cases: 

6 .  .  =  1  and  6 .  .  =  0 .  We  first  prove  (2.1)  for  the  case  6 .  .  =  1 .  Since  4>  is  increas- 
l j  l j  *  ii 

ing  in  its  arguments. 


sup{  t :  *(I(XU  >  t),  ....  I  (Xim  >  t))  =  1} 


(A4) 


<  sup{t:  *  (I(Xu>t) . I(X.m>t))  =  1). 

It  is  clear  that  the  left  side  of  (A4)  equals  and  the  right  side  of  (A4)  equals 

Y. . .  Hence 
ij 


T.  <  Y.  . . 
i  i) 


(A5) 


Since  6 .  .  =  1 , 
ij 


X.  .  <  T. . 
ij  i 


(A6) 


It  is  immediate  from  (A5)  and  (A6)  that  X„  <  Y  _  ,  which  implies  that  (2.1)  holds 

for  this  case.  We  now  prove  that  (2.1)  is  satisfied  if  =0.  Since  6 ^  =0,  it 

follows  that  X .  .  >  T.  =  Z .  .  .  Hence  0  =  4>.(I(X.1  >Z.  .),  ...,  I(X.  >  Z .  . ) )  .  Thus  it 
ij  l  ij  j  v  v  ii  ij  lm  1 2 

follows  from  (A2)  that 


Y.  .  <  Z.  .  . 
tl  i.l 


(A?) 


It  is  easv  to  see  that  (A5)  holds  for  this  case.  Thus  Y..=Z..,  which  implies 

ij  iJ 


that  (2.1)  is  satisfied  for  this  case. 
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